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ABSTRACT 


Observations  of  the  temporal  evolution  of  waves,  currents,  and  bed  response  data 
collected  by  an  instrumented  frame  deployed  in  12m  of  water  at  the  Monterey  Inner  Shelf 
Observatory  (MISO)  off  the  coast  of  Monterey,  California  are  analyzed  in  terms  of 
measured  wave  and  current  forcing  statistics  and  ripple  geometry.  During  the  year  2000, 
a  Broadband  Acoustic  Doppler  Current  Profiler  (BADCP)  collected  continuous  wave  and 
current  measurements.  Bed  morphology  was  continually  mapped  by  a  Scanning  Acoustic 
Altimeter  (SAA)  in  a  Im  alongshore  by  1.5m  cross- shore  area  immediately  offshore  from 
the  MISO  frame.  Relict  ripples  were  observed  to  dominate  the  bedforms  throughout 
much  of  the  year.  Ripple  growth  in  the  alongshore  direction  was  observed  during 
conditions  of  marginally  critical  flow  as  defined  by  the  critical  combined  wave  and 
current  Shields  parameter.  As  flow  conditions  increased  above  the  critical  level,  ripple 
growth  in  the  alongshore  direction  ceased,  and  cross- shore  wavelengths  began  to  grow 
and  dominate.  Together,  these  observations  and  data  sets  are  used  to  evaluate  the 
applicability  of  existing  ripple  prediction  algorithms.  Altogether,  five  models  are  tested, 
and  it  was  concluded  that  they  could  not  independently  predict  the  bed’s  response. 
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1.  INTRODUCTION 


Bottom  bedforms  develop  on  sandy  beds  in  shallow  coastal  waters  along  the  inner 
continental  self  in  response  to  flows  of  water  of  sufficient  strength  to  mobilize  the 
sediment.  As  these  bedforms  develop  and  ripples  are  created,  the  bottom  roughness 
changes,  altering  the  structure  of  the  turbulent  bottom  boundary  layer.  Thus,  ripple 
formation  has  a  significant  impact  on  bed  stress,  wave  dissipation  rates,  near  bed  velocity 
structure,  and  the  amount  of  sediment  suspended  in  the  water  column  [Grant  and 
Madsen,  1982].  To  better  understand  the  bottom  boundary  layer  hydrodynamics  and 
sediment  transport,  the  bottom  geometry  of  mobile,  sandy  beds  as  a  function  of  flow  and 
sediment  parameters  must  be  characterized. 

Numerous  laboratory  and  field  studies  have  been  conducted  to  determine 
relationships  between  wave  forcing,  grain  size,  and  ripple  geometry.  Most  of  these 
studies  have  focused  on  wave-formed  ripples  typical  of  confined  shallow  water 
environments  where  wave  stresses  dominate  over  those  caused  by  mean  currents. 
Laboratory  studies  pioneered  by  Bagnold  [1946]  and  continued  by  Yalin  and  Russell 
[1962],  Mogridge  and  Kamphuis  [1972],  and  Miller  and  Komar  [1980a],  among  others, 
have  sought  empirical  brmulae  to  classify  the  size  and  shape  of  sand  ripples  using  a 
variety  of  experimental  devices  including  flumes,  water  tunnels,  oscillatory  beds,  using 
different  sediment  sizes  and  densities. 

Bagnold  [1946]  originally  classified  wave -generated  ripples  into  two  categories: 
rolling  grain  ripples  and  vortex  ripples.  Rolling  grain  ripples  formed  as  sand  grains 
began  to  move  and  align  with  each  other  in  ridges  perpendicular  to  gentle  oscillatory 
flow.  Once  the  flow  velocity  increased  to  twice  that  of  the  onset  of  sediment  motion,  a 
vortex  formed  on  the  lee  side  of  the  ridge.  This  vortex  moved  sand  from  the  lee-side 
trough  to  the  crest,  effectively  increasing  the  ripple  steepness,  and  creating  vortex  ripples. 
Like  rolling  grain  ripples,  the  length  of  the  vortex  ripples  grew  until  equilibrium  was 
reached  with  the  wave  orbital  diameter,  or  washed  out  as  the  wave  stress  reached  sheet 
flow  conditions.  Bagnold  [1946]  concluded  that  vortex  ripples  were  the  primary  bed 
response  in  strongly  forced,  high  sediment  transport  flow  conditions. 
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Bagnold’s  1946  experiments  and  other  laboratory  studies  are  limited  by  the  level 
of  wave  conditions  generated.  Small-scale  laboratory  equipment  cannot  adequately 
reproduce  long  period,  high  energy  wave  conditions  often  experienced  in  the  field. 
Wave-dominated  field  studies  conducted  by  Dingier  [1974]  and  Miller  and  Komar 
[1980b]  took  into  account  a  wide  range  of  wave  heights  and  periods.  Field  studies  such 
as  these  are  complicated  when  the  wave  field  is  broadbanded  with  directional  spreading, 
making  relationships  between  wave  conditions  and  measured  ripple  geometry  more 
difficult  to  ascertain.  In  addition,  due  to  the  complexities  of  obtaining  concurrent 
measurements  of  wave  and  bed  characteristics,  field  studies  are  often  short  in  duration 
and  only  marginally  measure  the  bedforms. 

To  gain  a  more  complete  understanding  of  the  bed’s  response  to  wave  forcing, 
Nielsen  [1981]  combined  the  findings  of  laboratory  and  field  studies  to  model  ripple 
geometry  in  wave -dominated  environments.  He  developed  separate  expressions  for 
laboratory  and  field  conditions,  obtaining  formulae  for  ripple  wavelength  (X),  height  (rj), 
and  steepness  (pA).  For  laboratory  data,  he  obtained  the  ripple  wavelength  relationship 

Xl  a  =  2.2-0.345^  (1) 

where  a  is  the  water  semi-excursion  amplitude,  and  t|/w  is  the  wave  mobility  number 
defined  by 


=(aOi)^  l(s-l)gD 


(2) 


(0  is  the  radian  frequency  (27i/T,  where  T  is  the  wave  period),  s  is  the  relative  sediment 
density,  D  is  the  mean  grain  diameter,  and  g  is  the  gravity  acceleration  constant. 
Nielsen’s  [1981]  laboratory  derived  ripple  height  calculation  is 


r\la  =  0.215-0.022-J^ 


(3) 
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and  his  ripple  steepness  formula  is 


ri/X  =  0.182  -  0.240^^ 


(4) 


where  the  Shields  parameter  0  is  equal  to  one  half  the  mobility  number  multiplied  by 
Jonsson’s  [1967]  wave  friction  factor. 

For  field  data,  Nielsen  [1981]  concluded  that  the  ripple  wavelength,  height,  and 
steepness  for  quartz  sand  had  the  following  relationships 
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ri/X  =  0.342-0.34^ 


(7) 


respectively. 

In  their  evaluation  of  several  wave- dominated  field  studies,  Clifton  and  Dingier 
[1984]  identified  three  categories  of  ripples:  orbital,  anorbital,  and  suborbital.  Orbital 
ripples  form  when  the  ratio  of  wave  orbital  diameter  (dw)  to  grain  size  diameter  (D)  lies  in 
the  range  of  100-3000.  In  this  range.  Miller  and  Komar  [1980a]  concluded  that  orbital 
ripple  wavelength  (Xorb)  is  proportional  to  the  wave  orbital  diameter  according  to 


Krb  =  0.65d„ 


(8) 
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At  large  ratios  with  dw/D  (>5000),  ripple  spacing  becomes  nearly  independent  of 
wave  orbital  diameter.  These  ripples,  called  anorbital,  are  prevalent  in  fine  sandy  beds 
under  long  period  waves  and  have  ripple  wavelengths  which  scale  directly  with  a 
multiple  of  grain  size  [Clifton  and  Dingier,  1984]: 


X_==400D-600D 


(9) 


In  the  transitional  range  between  orbital  and  anorbital  ripples,  intermediate 
wavelength  features  called  suborbital  ripples  develop.  The  relationship  between  ripple 
wavelength  and  wave  orbital  diameter  /  grain  diameter  of  suborbital  ripples  varies. 

Wiberg  and  Harris  [1994]  analyzed  Clifton  and  Dingler’s  [1984]  ripple 
classification  system  for  both  laboratory  and  field  studies,  and  concluded  that  while  short 
wave  period  laboratory  equipment  produces  orbital  ripples,  ripple  conditions  on  the 
continental  shelf  most  often  fall  in  the  range  of  suborbital  or  anorbital.  By  analyzing 
several  sets  of  laboratory  and  field  ripple  measurements,  Wiberg  and  Harris  [1994] 
provide  different  estimates  of  orbital  and  anorbital  wavelength  according  to 


X 


orb 


0.62d„ 


(10) 


I 


ano 


=  535D 


(11) 


respectively.  The  average  ripple  steepness  for  orbital  and  anorbital  ripples  was 
determined  according  to  the  following  two  equations: 
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Wiberg  and  Harris  [1994]  contend  that  the  ratio  of  wave  orbital  diameter  to 
anorbital  ripple  height  best  classifies  ripples  into  orbital,  anorbital,  and  suborbital 
conditions  according  to  the  limits 
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While  equations  10-13  classify  ripple  wavelength  and  steepness  for  orbital  and 
anorbital  cases,  a  suborbital  wavelength  predictor  was  derived  using  the  boundaries 
established  in  equations  14(a)  and  14(c): 
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The  underlying  assumption  in  the  Nielsen  [1981]  and  Wiberg  and  Harris  [1994] 
models  is  that  the  oscillatory  flow  component  dominates  at  all  times.  This  is  not  always 
true,  as  flows  with  both  steady  and  oscillatory  components  are  commonly  present  in  the 
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continental  shelf  (most  noticeably  during  storms  ^lenn  and  Grant,  1987]),  and  can 
extend  to  the  seafloor  to  add  to  the  wave/oscillatory  motion  and  effectively  initiate 
sediment  motion  [Grant  and  Madsen,  1979].  Clifton  and  Dingier  [1984]  note  that  ripples 
become  asymmetric  in  the  direction  of  the  mean  current  under  oscillatory  flow  with  axial 
mean  currents. 

Until  recently,  there  have  been  few  investigators  of  bedforms  generated  from 
combined  wave  and  current  flow.  Laboratory  studies  conducted  by  Grant  and  Madsen 
[1979],  Glenn  and  Grant  [1987],  and  Khelifa  and  Ouellet  [2000]  as  well  as  field  studies 
by  Amos  et  al.  [1988]  and  Li  and  Amos  [1998]  all  present  empirical  formulae  for 
predicting  ripple  geometry  under  the  combined  actions  of  waves  and  currents.  Khelifa 
and  Ouellet  [2000]  conclude  that  empirical  formulae  used  to  model  ripple  wavelength 
and  height  for  pure  wave  motion  are  not  applicable  for  ripple  prediction  under  combined 
waves  and  currents.  Khelifa  and  Ouellet  [2000]  proceed  to  derive  empirical  formulae  for 
ripple  wavelength  and  height  from  laboratory  produced  ripple  measurements.  Their 
ripple  wavelength  equation  for  combined  flow  is 


—  =  1.9+0.081ii»(1+V.,)-0.741ii(1+>|/..) 


(16) 


where  the  combined  wave/current  orbital  diameter  designated  dew  is 
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with  0  representing  the  angle  between  waves  and  currents.  Mg  =  the  current  velocity  at  the 
top  of  the  boundary  layer,  and  T  =  the  wave  period,  and  the  combined  flow  mobility 
number 
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is  the  ratio  of  the  combined  flow  velocity  (JJcw)  squared  to  the  immersed  weight  per  unit 
area  of  the  bed  according  to 
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with  the  depth- averaged  current  velocity  m  replacing  Wg  from  (17). 

Khelifa  and  Ouellet’s  [2000]  ripple  height  formula  is  expressed  as 


^  =  0.32  +  0.017h^(l+V,J-0.142h(l+v,J  (20) 


Similarly,  Li  and  Amos  [1998]  applied  wave-dominated  ripple  prediction 
equations  to  their  field  data  and  concluded  that  wave-ripple  predictors  significantly  over¬ 
predict  ripple  height  and  roughness  for  combined  flows.  Li  and  Amos  [1998]  classified 
combined- flow  ripples  into  three  categories  based  on  the  ratio  of  skin  friction  wave  shear 
{uws)  to  skin  friction  current  shear  velocity  (ucs)'-  wave  dominant  (Uws/Ucs  >  1-25),  current 
dominant  ^w/Ucs  <  0.75),  and  combined  wave-current  ripples  (0.75  <  Uw/Ucs  <  1.25). 
Their  proposed  ripple  height  predictors  for  both  wave  dominant  and  combined  flow 
produced  ripples  are  as  follows. 

For  wave  dominant  ripples. 
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where  Ucwb  is  the  bedload  shear  velocity  and  Ucr  is  the  critical  shear  velocity  for  bedload 
transport.  The  wave- current  ripples  were  modeled  according  to 
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In  addition  to  waves  and  mean  currents,  other  oceanic  phenomena  such  as  bound 
long  waves  affect  sediment  transport.  Bound  long  waves  are  within  the  infragravity 
frequency  range  (f  <  0.05  Hz,  or  T  >  20s),  and  serve  as  a  restoring  force  for  the  variation 
in  short  wave  momentum  of  incident  wave  groups.  As  a  wave  train  consisting  of  wave 
groups  (surf  beat)  and  the  associated  bound  long  waves  begin  to  shoal,  the  waves  become 
non-linear  with  peaked  asymmetric  crests  and  long,  flat  troughs  [Longuet -Higgins  and 
Stewart,  1962].  The  waves’  positively  skewed  orbital  velocity  suspends  and  transports 
sediment  onshore,  while  the  lower  frequency  currents  transport  suspended  sediment 
either  onshore,  offshore,  of  alongshore  depending  on  the  current’s  mean  direction.  The 
bound  long  waves  under  wave  groups,  however,  provide  a  net  offshore  transport.  Both 
Shi  and  Larsen  [1984]  and  Deigaard  et  al.  [1999]  conclude  that  because  bound  long 
waves  produce  an  offshore  flow  under  high  waves  in  the  group  when  the  suspended 
sediment  concentrations  are  high,  and  an  onshore  flow  under  the  low  waves  with  low 
suspended  sediment  concentrations,  a  net  offshore  sediment  transport  results.  An 
investigation  of  the  relative  importance  of  these  different  transport  mechanisms  derived 
from  laboratory  and  field  studies  would  greatly  aid  in  predicting  net  sediment  transport 
[Deigaard  et  al.,  1999]. 

To  gain  a  better  understanding  of  sediment  transport  and  ripple  geometry,  several 
methods  of  observing  ripple  structure  and  migration  have  been  utilized  in  the  field. 
Many  previous  field  studies  have  relied  upon  crude  in  situ  observations  conducted  by 
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divers,  which  greatly  limits  temporal  and  spatial  sampling  of  bottom  bedforms.  Boyd  et 
al.  [1988],  Amos  et  al.  [1988],  and  Wheatcroft  [1994]  used  photographic  measurement 
time- series  to  characterize  ripple  wavelength,  orientation,  and  migration  rates  under 
varying  oceanographic  conditions.  Li  and  Amos  [1998]  also  used  photographic 
equipment  to  monitor  the  bed,  but  included  a  shadow  bar  to  aid  in  determining  ripple 
height.  These  optical  methods  for  two-dimensional  bedform  monitoring  lack  accuracy, 
and  perform  poorly  during  periods  of  high  suspended  sediment  concentration. 

Hay  and  Wilson  [1994]  and  Traykovski  et  al.  [1999]  employed  rotary  sidescan 
sonar  systems  to  quantify  bedform  evolution  in  their  nearshore  field  studies.  High 
resolution  images  of  bedform  development  were  obtained  in  this  manner,  and  ripple 
wavelength  and  migration  rates  were  determined.  Similar  to  the  photographic  systems, 
however,  this  apparatus  is  limited  by  its  inability  to  ascertain  bedform  amplitudes. 

The  Monterey  Inner  Shelf  Observatory  (MISO)  was  designed  and  implemented 
by  Stanton  [1999]  to  obtain  continuous  observations  of  bedforms  in  12m  of  water  while 
simultaneously  monitoring  wave  and  current  conditions.  A  scanning  acoustic  altimeter 
and  an  acoustic  Doppler  current  profiler  obtained  high  resolution  measurements  of  ripple 
geometry  and  flow,  respectively.  The  objective  of  this  paper  is  to  analyze  the  wave, 
current,  and  bedform  response  data  collected  by  MISO  off  the  coast  of  Monterey, 
California  during  the  year  2000.  Due  to  the  length  of  the  time- series,  the  data  spans  a 
wide  range  of  dynamic  oceanic  conditions  and  various  bed  responses,  spanning  relict, 
active  ripple  formation,  and  sheet  flow  conditions.  Existing  ripple  models  are  evaluated 
using  these  field  data. 
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II.  SITE  DESCRIPTION 


The  Monterey  Inner  Shelf  Observatory  (MISO)  is  part  of  an  on-going  program 
focused  on  littoral  oceanographic  observation  and  modeling  at  the  Naval  Fbstgraduate 
School,  Monterey,  CA.  MISO,  a  component  of  the  Rapid  Environmental  Assessment 
Laboratory  (REAL),  consists  of  an  underwater  frame  supporting  five  instruments.  It  is 
located  in  12m  of  water  600  m  off  the  Del  Monte  Beach  at  the  southern  end  of  the 
Monterey  Bay  (Eat.  36°  36.29'  N;  Long.  121°  052.53'  W;  Eigure  1)  on  a  broad,  sandy 
channel  between  shale  beds  to  the  north  and  south.  The  frame  is  mounted  over  a  mobile 
bed  of  well- sorted  quartz  sand  with  a  0.260  mm  median  grain  diameter  (dso)  (Eigure  2). 
Waves  incident  on  the  MISO  location  tend  to  be  narroxwbanded  from  the  northwest  due 
to  the  strong  refractive  effects  of  the  Monterey  Bay  Canyon  as  well  as  the  protection 
provided  by  the  Point  Pinos  and  Santa  Cruz  headlands  (Eigure  1). 

MISO  is  connected  to  a  shore  terminus  by  a  fiber  optic  data  and  power  cable.  A 
meteorological  station  in  the  sand  dunes  inshore  augments  the  data  collected  underwater 
by  measuring  10m  height  wind  speed  and  direction,  air  temperature  and  dew  point,  short 
and  long  wave  incident  radiation,  barometric  pressure,  and  rain  fall  rates.  Eigure  3 
illustrates  the  MISO  and  meteorological  station’s  configuration. 

MISO  began  collecting  data  on  July  26,  1999,  and  the  meteorological  station 
began  on  April  28,  2000.  Although  the  data  set  is  not  complete  (due  to  routine 
maintenance  and  instrument  removal/repair),  a  large  portion  of  continuously  sampled 
data  is  available  for  study.  This  study  focuses  exclusively  on  the  year  2000  data  as  it 
contains  the  most  continuous  time- series  of  velocity  and  bedform  measurements  (Eigure 
4).  This  time  period  provides  comprehensive  data  spanning  energetic  winter  to  relatively 
quiescent  summer  conditions.  The  average  significant  wave  height  for  this  time- series  is 
0.5  m  with  a  peak  around  2  m. 
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Location  of  the  MISO  instrument  platform  off  Monterey,  CA.  The  star  represents 
the  location  of  the  MISO  array. 


MISO  Sand  Gran  Size 


Figure  2.  Sand  grain  size  distribution  of  six  sand  samples  taken  at  the  MISO  site.  The 
median  grain  size  (dso)  of  these  samples  is  0.260  mm. 
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MISO  Instrument  Timelines 
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Figure  4.  MISO  instrument  timeline  for  the  year  2000.  The  above  timeline  indicates  when 
each  instrument  mounted  on  the  MISO  frame  and  meteorological  tower  collected 
data.  Gaps  indicate  instrument  maintenance  periods.  [MET  -  meteorological 
station  instruments;  SURFCAM  -  surf  camera  mounted  on  the  meteorological 
station;  BCDVl  -  Bistatic  Coherent  Doppler  Velocity  and  Sediment  Profiler; 
SANDCAM  -  Structured  Light  and  Camera;  XYALT  -  Scanning  Acoustic 
Altimeter  (SAA);  PARO  -  Paro- Scientific  pressure  sensor;  BADCP  -  Bottom 
Acoustic  Doppler  Current  Profiler] 
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III.  INSTRUMENTATION  AND  DATA  PROCESSING 


The  MISO  frame  mounts  five  instruments  for  monitoring  bedforms  and  wave  and 
current  forcing:  a  Broadband  Acoustic  Doppler  Current  Profiler  (BADCP),  Paro- 
scientific  pressure  sensor,  Scanning  Acoustic  Altimeter  (SAA),  Structured  Light  and 
Camera  (SLAC),  and  a  Bistatic  Coherent  Doppler  Velocity  and  Sediment  Profiler 
(BCDVSP).  Figure  5  illustrates  the  placement  of  these  instruments  on  the  MISO  frame. 
The  BADCP  measures  the  current  velocities  every  Im  from  the  frame  to  the  near  ocean 
surface,  allowing  mean  currents  and  characteristics  of  surface  waves  to  be  determined. 
The  Paro- scientific  pressure  sensor  measures  the  height  of  water  above  the  frame  to 
produce  a  record  of  tides  and  long  period  waves.  The  SAA  continuously  maps  the  sandy 
bed  offshore  of  the  MISO  frame,  producing  a  time- series  of  the  changing  bedform 
conditions  for  a  Im  alongshore  by  1.5m  cross-shore  area,  while  the  SLAC  produces 
photographic  images  of  a  Im  by  Im  area  to  quantify  small-scale  bedforms.  The 
BCDVSP  measures  the  velocity  and  sediment  concentration  above  the  bed.  These  small- 
scale  measurements  are  used  to  define  the  mean  and  turbulent  structure  of  the  water 
column  above  the  bed  resulting  from  combined  flow  conditions  over  the  evolving 
bedforms.  This  study  utilizes  data  from  the  BADCP  and  SAA  as  they  provide  the  most 
continuous  velocity  and  ripple  measurements,  respectively  (Figure  4),  for  a  year-long 
period. 


A.  BROADBAND  ACOUSTIC  DOPPLER  CURRENT  PROFILER  (BADCP) 

I.  Data  Acquisition 

An  RDI  Instruments’  Broadband  Acoustic  Doppler  Current  Profiler  (BADCP) 
measured  horizontal  current  velocities  at  the  MISO  site  by  transmitting  broadband  pulses 
of  1.2  MHz  frequency  sound  waves  along  each  of  four  narrow  acoustic  beams,  slanted  20 
degrees  from  the  vertical.  Each  transducer  then  measured  the  echoes  returned  from 
sound  scatterers  (small  particles  or  plankton)  suspended  in  the  water  column  at  Im  range 
intervals  up  to  the  sea  surface.  As  the  mean  currents  carried  the  sound  scatterers  to  or 
from  the  ADCP,  the  sound  waves  at  the  scatterers’  location  become  Doppler- shifted  to 
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STRUCTURED  LIGHT 


instruments  include  (from  left  to  right):  camera  component  of  the  Structured  Light 
and  Camera  (SLAC),  Bistatic  Coherent  Doppler  Velocity  and  Sediment  Profiler 
(BCDVSP),  structured  light  component  of  the  SLAC,  Scanning  Acoustic 
Altimeter  (SAA),  Paro- scientific  pressure  sensor,  and  the  Broadband  Acoustic 
Doppler  Current  Profiler  (BADCP). 


higher  or  lower  frequencies.  This  Doppler- shift  is  proportional  to  the  relative  velocity 
between  the  ADCP  and  the  scattering  particles.  By  measuring  Doppler  shifts  with  a  four- 
beam  configuration,  the  ADCP  resolved  horizontal  current  velocities  into  Im  range - 
averaged  depth  cells,  or  bins,  at  a  sampling  frequency  of  0.976  Hz.  A  high  resolution 
strain  gage  pressure  transducer  in  the  BADCP  head  measured  the  water  pressure  at  the 
same  sample  frequency. 

The  horizontal  ^elocities  measured  at  the  MISO  site  were  oriented  relative  to  the 
ADCP  instrument  and  MISO  frame  at  a  magnetic  heading  of  313°  M.  Following  an 
error-removal  process,  the  current  velocity  measurements  were  remapped  into  a  cross¬ 
shore  normal  coordinate  system.  The  velocity  vectors  were  rotated  7°  -i-  90°  such  that  the 
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Figure  6.  Cross- shore  normal  coordinate  system  for  horizontal  current  velocities  measured 

by  the  Bottom  Acoustic  Doppler  Current  Profiler  (ADCP)  as  well  as  the 
orientation  for  the  Scanned  Acoustic  Altimeter  (SAA)  at  the  MISO  site  off  Del 
Monte  Beach,  Monterey,  CA. 
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+u  vector  corresponded  to  the  nominal  cross-shore  (CS)  direction  and  the  -tv  vector 
corresponded  to  the  alongshore  (AS)  direction  as  shown  in  Figure  6. 

The  raw  ADCP  data  time  series  was  segmented  into  one -hour  ensembles  for 
statistical  analysis  to  provide  stable  estimates  of  dynamic  forcing  variables  at  the  site. 
The  first  Im  bin  of  the  ADCP  vertical  profile  was  discarded  due  to  velocity  errors 
induced  by  the  water’s  interaction  with  MISO’s  equipment  frame.  The  second  velocity 
bin  (~2m  off  the  bottom,  or  z=10m)  was  utilized  to  characterize  the  near-bed  velocity 
field  well  outside  of  the  thin,  oscillatory  boundary  layer.  The  linear  wave  transfer 
functions  for  swell  waves  typical  of  the  site  are  presented  in  Figure  7.  The  near-bed 
depths  of  12m  produced  transfer  function  values  within  5%  of  those  at  10m  depth  for 
typical  wave  periods  of  8,  10,  and  15s.  For  this  reason,  the  10m  velocity  measurements 
were  used  to  represent  the  near-bed  velocities  just  outside  the  thin  oscillatory  boundary 
layer  above  the  bed. 

2.  Error  Removal  Procedure 

A  concern  in  processing  the  multiyear  ADCP  time- series  from  the  MISO 
experiment  was  distinguishing  the  true  current  velociies  near  the  bed  from  occasional 
erroneous  values.  Even  infrequent  spurious  velocity  estimates  had  to  be  removed  to 
ensure  meaningful  high- order  statistics.  Errors  not  representative  of  true  current 
velocities  likely  occurred  in  the  broadband  ADCP  data  as  a  result  of  fish  or  seaweed 
passing  in  front  of  the  transducers,  which  in  turn  generated  poor  beam  correlation  values. 
These  intermittent  errors  could  be  seen  upon  visual  inspection  of  the  raw  time- series  plots 
of  alongshore  (AS)  and  cross-shore  (CS)  current  velocity  measurements.  Error  detection 
and  removal  using  an  interactive  editor,  however,  would  have  been  tedious  and  time 
consuming.  Instead,  these  erroneous  data  points  were  identified  and  removed  in  an 
automated  error  removal  routine  utilizing  four  discriminating  steps:  (1)  a  correlation  level 
test,  (2)  a  normalized  difference  test,  (3)  a  temporal  difference  test,  and  (4)  an  outlier 
value  test.  All  current  velocity  data  were  deglitched  in  hourly  sections  prior  to  further 
analysis  using  the  following  criteria. 
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Figure  7.  Vertical  profile  of  linear  wave  theory  transfer  function  for  ADCP  velocity 

measurements  of  5,  8,  10,  and  15s  period  waves  in  12m  of  water.  The  transfer 

function  ^ ,  where  k  is  the  wavenumber  and  h+z  is  the 

cosh(A:(/f+  z)) 

distance  above  the  bed,  produced  roughly  equivalent  values  for  10m  and  12m 
depth  velocity  measurements. 
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a. 


Step  1:  ADCP  Correlation  Level  Test 


In  addition  to  ascertaining  the  along-beam  current  velocity  at  a  given  time 
and  given  range  bin,  the  ADCP  also  recorded  the  data  quality  returned  in  the  form  of  a 
correlation  value.  The  correlation  magnitude  measured  the  normalized  echo 
autocorrelation  at  the  lag  used  for  estimating  the  Doppler  phase  change  of  each  beam 
(RDI,  1998).  The  Workhorse  ADCP  represented  correlation  values  in  a  linear  scale 
between  0  and  255,  with  255  being  a  perfect  correlation  (i.e.  a  solid  target  with  no  system 
or  acoustic  noise),  and  0  indicating  the  poorest  correlation.  Low  correlation  values  in  the 
MISO  ADCP  data  may  have  been  caused  by  low  scattering  particle  concentrations  or 
near-field  objects  obscuring  the  transducer,  resulting  in  poor  signal- to- noise  ratios.  To 
remove  erroneous  velocity  estimates  resulting  from  poor  correlation  magnitudes,  all  four 
beams’  correlation  values  were  first  added  together  at  each  time  and  range  bin.  Upon 
determining  the  mean  and  standard  deviation  of  the  summed  correlation  values  over  a 
one- hour  interval,  a  threshold  was  established: 


C,Hr=C^,an  -  CoVrC  *  C 


(23) 


where  C  represents  the  correlation  time- series,  and  CorrC  is  a  constant  pre-selected  value 
(see  Table  1).  If  the  correlation  values  dropped  below  this  threshold  at  any  point  in  the 
one-hour  time- series,  the  velocity  data  was  flagged  as  erroneous,  and  removed  from  the 
data  set.  This  method  provided  adaptive  criteria  for  rejecting  velocity  estimates  using  the 
correlation  time- series,  which  could  vary  slowly  as  backscatter  and  fluid  flow  conditions 
changed  through  the  yearlong  time- series. 
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b.  Step  2:  Normalized  Range  Bin  Difference  Test 

A  clear  indication  of  errors  or  spikes  in  the  data  was  observed  when  the 
time-series  of  two  adjacent  velocity  bins  were  compared.  While  it  was  expected  that  the 
velocity  measurements  between  each  bin  would  vary  at  any  given  time,  the  along-beam 
difference  velocities  are  bounded  for  a  given  wave  height  above  the  bed  and  mean  current 
shear.  Large  differences  in  velocity  from  one  bin  to  another  near  the  bed  might  indicate 
the  presence  of  fish  or  other  anomalous  scattering  particles.  In  order  to  detect  significant 
vertical  differences  of  horizontal  current  velocity,  the  difference  between  two  bins’ 
velocity  vectors  was  normalized  by  their  standard  deviation.  If  the  normalized  difference 
between  the  bins  at  each  sample  point  was  greater  than  NormC,  the  velocity  data  was 
rejected: 


U 


thr 


U 


difference 


U 


rms 


>  NormC 


(24) 


where  Udijjerence  =  Ubin2-Ubin3,  U rms  Is  Calculated  over  each  one- hour  interval,  and  NormC 
is  a  predetermined  threshold  value,  provided  a  conservative  threshold  that  adapted  to 
the  average  wave  forcing  over  each  hour. 


c.  Step  3:  Temporal  Difference  Test 

The  third  procedure  in  the  deglitching  sequence  removed  errors/spikes  not 
identified  in  the  correlation  or  normalized  difference  correction  methods  by  testing  the 
fluid  acceleration  between  successive  samples.  First  differences  in  time  were  calculated 
for  the  velocity  time- series,  and  values  exceeding  the  threshold  ±DespC  *  a  and  not 
extending  for  more  than  4  sample  intervals  (~4  seconds)  were  replaced  with  interpolated 
values  using  a  linear  interpolation  scheme.  The  objective  of  this  criterion  was  to  remove 
large  slew,  single  point  outliers  from  the  time- series. 
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d.  Step  4:  Outlier  Value  Test 

As  a  final  precaution,  all  remaining  outlying  data  points  were  removed  by 
first  demeaning  the  velocity  data  set.  All  demeaned  velocity  data  points  whose  absolute 
value  was  greater  than  the  large  pre-selected  threshold  value  OutC  *  o  were  removed  and 
replaced  with  flagged  error  values: 


Othr=\Udemeaned\>OutC*U 


(25) 


Once  all  of  the  errors  were  identified  and  replaced  with  flagged  values,  a 
gap  replacement  routine  filled  gaps  <  5  time  intervals  (~5  seconds)  in  length  with  the 
mean  current  velocity.  Thus,  small  gaps  resulting  from  errors  induced  by  low  correlation 
values  or  spurious  scattering  particles  in  the  water  column  would  not  affect  subsequent 
statistical  measurements.  Larger  gaps,  however,  remained  after  the  replacement  process 
and  posed  a  problem  in  higher-order  statistical  computations.  For  this  reason,  high-order 
statistical  quantities  were  only  derived  from  continuous  data  segments  greater  than 
twenty  minutes  in  length. 


e.  Statistical  Analysis  of  the  Error  Removal  Procedure 

Great  care  was  taken  to  make  the  data  set  as  self-consistent  as  possible  by 
establishing  conservative  threshold  levels  that  would  not  reject  valid,  large  amplitude 
data  points  in  the  velocity  time- series.  The  alongshore  current  (AS)  velocity  time- series 
of  yeardays  2000078  (MAR  2000)  and  2000137  (MAY  2000)  are  used  as  examples  to 
illustrate  the  removal  of  erroneous  velocity  samples  and  the  choice  of  threshold  values. 
Figure  8  shows  that  day  078  was  highly  energetic  with  AS  velocities  2m  above  the  bed 
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Ycarday  2000078:  Error  Removal  (AS  Component) 
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Figure  8.  Alongshore  (AS)  raw  and  processed  ADCP  velocity  data  time- series  for  high- 
energy  yearday  2000078. 

oscillating  between  -0.2  and  0.2  m/s,  while  Figure  9  shows  that  day  137  was  calmer  as 
evidenced  by  lower  velocity  magnitudes  ranging  from  -0.1  to  0.1  m/s  at  the  same  depth. 
The  corresponding  cross-shore  component  magnitudes  (not  pictured)  ranged  from  -0.8  to 
0.8  m/s  and  -0.2  to  0.2  m/s  for  yeardays  078  and  137,  respectively.  The  processed  data 
in  both  figures  is  a  result  of  applying  conservative  threshold  levels  to  the  raw  ADCP 
velocity  data  at  each  stage  of  the  error  removal  process.  Table  1  provides  a  summary  of 
the  various  threshold  categories  tested  for  each  step  of  the  error  removal  process,  and 
indicates  the  selected  threshold  levels  used  to  obtain  the  processed  data  in  these  figures. 

To  measure  the  effectiveness  of  the  error  removal  sequence,  kurtosis 
coefficients  were  calculated  for  various  threshold  levels  at  each  hour  time  interval.  The 
kurtosis  coefficient,  or  4'*’  central  moment  divided  by  the  fourth  power  of  the  standard 
deviation,  is  a  measure  of  distribution  peakedness,  and  is  highly  sensitive  to  erroneous 
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Ycarday  2000137:  Error  Removal  (AS  Component) 


-02 


■04 


137,1 


I  :  I  I  l  I  I 

137  2  137  3  1374  137  3  137  6  137,7  137  0  137  9 

Yearday  -  2000 


Figure  9.  Alongshore  (AS)  raw  and  processed  ADCP  velocity  data  time- series  for  low- 
energy  yearday  2000137. 

Table  1.  Error  removal  threshold  names  and  values.  Selected  threshold  values  correspond 
to  the  levels  used  in  the  error  removal  process,  while  low  and  high  values  are  used 
to  validate  those  selections. 


Threshold 
Correction  Name 

Threshold 

Abbreviation 

Low 

Threshold 

Values 

Selected 

Threshold 

Values 

High 

Threshold 

Values 

Correlation 

CorrC 

3 

4 

5 

Normalized 

Range -Bin 
Difference 

NormC 

4 

5 

6 

Temporal 

Difference 

Despike 

DespC 

3 

4 

5 

Outlier 

OutC 

3 

4 

5 
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Y8B«Ja/  21X>0078:  AS  Current  Velocity  Kurtosis 
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Figure  10.  Hourly  kurtosis  coefficients  for  alongshore  (AS)  ADCP  velocity  data  time-series 
for  yearday  2000078:  (a)  using  low  threshold  values,  (b)  using  the  selected 
threshold  values,  and  (c)  using  high  threshold  values. 


outlying  data  points.  For  this  reason,  it  was  used  to  determine  the  robustness  of  the  error 
removal  process.  A  kurtosis  value  of  3  would  indicate  that  the  velocity  distribution  was 
normal,  while  values  greater  than  3  would  indicate  outlying  data  points.  If  the  kurtosis 
values  dropped  below  3,  the  data  distribution  became  flatter  than  normal,  indicating  too 
many  data  points  had  been  removed  from  the  time- series.  The  automated  data  editing 
algorithm  required  that  true  changes  in  high  order  statistics  be  retained  while  occasional 
erroneous  values  are  removed.  Figures  10  and  11  show  the  AS  hourly  kurtosis  values  for 
each  example  yearday  with  different  threshold  categories:  low,  selected,  and  high. 

Thresholds  were  systematically  chosen  for  each  deglitching  criterion  by 

observing  the  number  of  data  points  being  removed  or  replaced  as  well  as  the  data 

distribution’s  peakedness  response  to  various  correction  levels.  The  threshold  levels  that 

were  conservatively  selected  are  as  follows:  correlation  correction  (CorrC)  =  4, 
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Y«B«Ja/  21)00137:  AS  Current  Velocity  Kurtosis 


Figure  11.  Hourly  kurtosis  coefficients  for  alongshore  (AS)  ADCP  velocity  data  time- series 
for  yearday  2000137:  (a)  using  low  threshold  values,  (b)  using  the  selected 
threshold  values,  and  (c)  using  high  threshold  values. 


normalized  range  bin  difference  correction  (^ormC)  =  5,  temporal  difference  despike 
correction  (DespC)  =  4,  and  outlier  correction  (OutC)  =  4.  Figures  12  and  13  illustrate 
how  many  data  points  were  removed  /  replaced  according  to  each  step  with  the  selected 
threshold  levels  for  both  yearday s.  The  number  of  corrected  data  points  for  each  criterion 
as  well  as  the  total  number  of  data  points  is  contained  in  Table  2.  The  results  of  lowering 
the  threshold  values  (i.e.  removing  more  data  points  by  decreasing  each  selected  level  by 
1)  as  well  as  the  number  of  data  points  removed  /  replaced  in  accordance  to  high 
threshold  levels  (i.e.  retaining  more  data  points  by  increasing  each  selected  level  by  1)  are 
also  displayed  in  Table  2  and  plotted  in  Figures  12  and  13. 
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Vearday  2000078:  Error  Removal  Proceduro  (AS  Component) 
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Figure  12.  Data  points  flagged  as  errors  for  the  alongshore  (AS)  ADCP  velocity  data  on 
yearday  2000078  according  to  each  step  of  the  error  removal  procedure: 

(a)  correlation  correction,  (b)  normalized  range  bin  difference  correction, 

(c)  temporal  difference  despike  correction,  and  (d)  outlier  correction  for  low, 
selected,  and  high  threshold  values. 

One  consequence  of  lowering  each  threshold  value  to  CorrC  =  3,  NormC 
=  4,  DespC  =  3,  and  OutC  =  3  was  a  significant  reduction  in  kurtosis  values  due  to  the 
removal  of  too  many  data  points  (Figures  10(a)  and  11(a)).  For  energetic  days,  reducing 
the  correlation  correction  from  4  to  3  increased  the  number  of  data  points  removed  by  a 
factor  of  5-10,  while  reducing  the  temporal  difference  and  outlier  correction  levels 
produced  similar  results.  For  low-energy  days,  a  reduction  in  the  correlation  correction 
threshold  was  not  as  significant.  Reducing  the  temporal  difference  and  outlier  correction 
thresholds,  however,  increased  the  number  of  data  points  rejected  ten- fold  in  low-energy 
days.  Figures  12  and  13  both  show  too  many  data  points  being  identified  as  erroneous  in 
the  low  threshold  time- series.  This  conclusion  is  validated  upon  inspection  of  the 
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Vearday  2000137:  Eror  Removal  Procedure  (AS  Cornponent) 
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Figure  13.  Data  points  flagged  as  errors  for  the  alongshore  (AS)  ADCP  velocity  data  on 
yearday  2000137  according  to  each  step  of  the  error  removal  procedure: 

(a)  correlation  correction,  (b)  normalized  range  bin  difference  correction, 

(c)  temporal  difference  despike  correction,  and  (d)  outlier  correction  for  low, 
selected,  and  high  threshold  values. 


corresponding  kurtosis  calculations.  All  low  threshold  values  produced  current  velocity 
kurtosis  values  lower  than  the  normal  distribution  value  of  3.  Clearly,  removing  too 
many  data  points  resulted  in  a  significantly  altered  velocity  distribution. 

Increasing  the  threshold  values  to  CorrC  =  5,  NormC  =  6,  DespC  =  5,  and 
OutC  =  5,  on  the  other  hand,  proved  to  retain  too  many  erroneous  outlying  data  points 
(particularly  during  the  low-energy  days  such  as  day  137).  While  the  goal  in  removing 
errors  from  the  data  set  was  to  retain  as  much  of  the  large  amplitude  data  points  as 
possible,  these  high  thresholds  were  retaining  many  outlying  data  points  and  producing 
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Table  2.  Number  of  erroneous  data  points  for  the  alongshore  ADCP  velocity  data  on 
yearday  2000078  and  2000137  according  to  each  step  of  the  error  removal 
procedure  using  different  threshold  level  categories. 


Day  078 

Low 

Selected 

High 

Total  Nuber  of 

CorrC 

257 

53 

8 

Data  Points 

NormC 

19 

3 

2 

67077 

DespC 

0 

0 

0 

OutC 

173 

5 

0 

Total 

449 

61 

10 

Day  137 

Low 

Selected 

High 

Total  Nuber  of 

CorrC 

141 

16 

3 

Data  Points 

NormC 

36 

21 

12 

67311 

DespC 

486 

44 

3 

OutC 

218 

33 

13 

Total 

881 

114 
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data  distributions  with  high  kurtosis  values  close  to  those  of  the  raw  data  distributions 
(Figures  10(c)  and  11(c)). 

For  these  reasons,  the  selected  threshold  values  were  established  at  CorrC 
=  4,  NormC  =  5,  DespC  =  4,  and  OutC  =  4.  Of  particular  interest  was  the  fact  that  the 
high-energy  days  contained  more  correlation  corrections  than  the  low-energy  days.  The 
numbers  of  normalized- difference,  temporal  difference,  and  outlier  corrections,  however, 
were  more  frequent  during  low-energy  days.  The  high-energy  days  contained  raw 
velocity  data  kurtosis  values  close  to  3,  while  the  low-energy  days  had  very  large,  erratic 
kurtosis  peaks  resulting  from  data  errors  in  the  raw  data  distribution.  Both  of  these 
observations  are  reasonable.  Very  energetic  days  with  swift  oscillatory  currents 
generated  low  correlation  values  due  to  high  fluid  speeds  through  the  sample  volumes, 
and  contained  relatively  few  outlying  data  points.  (Fish  are  less  likely  to  be  grazing  near 
the  bed  during  high  currents.)  In  contrast,  low-energy  days  with  weak  current  velocities 
reduced  the  possibility  for  low  correlation  values  by  making  occasional  detections  of  fish 
more  noticeable.  Any  scattering  particles  in  the  water  column  had  a  higher  residence 
time  in  the  ADCP’s  beam  path  than  during  very  energetic  conditions.  This  would 
account  for  the  increase  in  outlying  data  points  and  the  consequent  peaks  in  kurtosis 
values  obtained  during  low-energy  time  periods. 
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To  better  demonstrate  the  effectiveness  of  the  selected  thresholds,  Figure 
14  presents  the  error-removal  sequence  for  one  hour  of  yearday  2000176  velocity  data 
where  there  were  large  spikes  in  kurtosis  values  for  the  raw,  unprocessed  data.  The  need 
for  these  four  largely  independent  criteria  is  also  illustrated.  On  yearday  2000176,  hour 
13,  the  kurtosis  coefficient  for  the  AS  velocity  component  was  ~  16.7.  Figure  14(a) 
shows  the  time- series  for  this  hour.  Three  peaks  (one  negative,  two  positive)  of  outlying 
data  points  in  the  raw  data  clearly  show  why  the  kurtosis  coefficient  was  high.  Both  the 
correlation  correction  (Figure  14(b))  and  the  normalized  range  bin  difference  correction 
(Figure  14(c))  removed  two  errors  using  the  selected  threshold  values.  The  despiking 
routine  and  the  outlier  correction  both  resolved  the  outliers  near  176.565  by  removing 
and  interpolating  across  the  spike  using  a  selected  threshold  value  of  4o  (Figure  14(d)). 
The  resulting  kurtosis  coefficient  from  the  error  removal  sequence  is  ~  3.00.  The 
deglitching  sequence  returned  the  AS  velocity  distribution  of  that  hour  to  normal  while 
removing  and  replacing  only  29  /  2796  data  points  (1.04%).  If  that  same  hour  of  data 
was  deglitched  according  to  low  threshold  values,  50  /  2796  (1.79%)  of  the  data  would  be 
rejected  and  corrected  producing  a  kurtosis  value  of  2.88.  High  threshold  standards 
would  remove  only  21  /  2796  (1.07%)  of  the  data,  yielding  a  kurtosis  coefficient  of  3.13. 
The  selected  threshold  values  produced  the  best  data  distribution  with  minimal 
corrections  for  the  yearlong  time- series  for  twelve  audit  cases  chosen  from  the  yearlong 
time- series. 


3.  Statistical  Analysis 

Forcing  of  the  bed  was  characterized  by  calculating  statistics  of  the  2m  level 
current  time- series  for  each  one -hour  interval  through  the  year  2000.  The  mean  and  the 
RMS  velocities  were  the  only  quantities  calculated  on  the  raw  data  time- series.  The 
RMS  velocity  time- series  provided  a  good  indication  of  energetic  wave  forcing  yeardays. 
The  higher  order  statistics  skewness  and  kurtosis  were  calculated  from  the  error- removed 
gapless  time- series.  The  skewness  coefficient  is  a  measure  of  the  degree  of  asymmetry  of 
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Figure  14.  Error-removal  sequence  for  the  thirteenth  hour  of  yearday  2000176  ADCP 
velocity  data:  (a)  raw  and  processed  data  using  the  selected  threshold  levels, 

(b)  ADCP  beam  correlation  values  with  a  4o  threshold  level,  (c)  normalized  range 
bin  difference  using  a  cutoff  value  of  5,  and  (d)  temporal  difference 
de spike/outlier  corrected  data  using  a  4a  threshold  level. 
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the  velocity  distribution,  and  is  calculated  by  dividing  the  third  moment  by  the  cube  of 
the  standard  deviation.  A  skewness  coefficient  of  0  indicated  equal  oscillatory  flow  in 
either  direction.  A  negative  skewness  (in  the  cross- shore  direction)  indicates  a  net 
offshore  flow,  while  a  positive  skewness  coefficient  represents  net  on-shore  flow.  The 
kurtosis  coefficient,  or  the  ratio  of  the  fourth  moment  to  the  variance  squared,  represents 
the  peakedness  or  flatness  of  the  velocity  component  time- series. 

Using  the  time- series  provided  by  the  pressure  sensor  in  the  BADCP  and  by 
applying  a  linear  wave  theory  transfer  function  to  the  surface,  significant  wave  height 
{Hi/s)  was  calculated  for  each  hour  using  the  following  equation: 


H 


1/3 


=  4*0 


pressure 


(26) 


where  ^pressure  is  the  standard  deviation  of  the  pressure  data  for  one  hour  scaled  to  m  of 
water.  Each  hour’s  averaged  pressure  data  was  also  used  to  monitor  tidal  fluctuations  at 
the  MISO  site. 

A  Maximum  Entropy  Method  (MEM)  directional  wave  spectra  algorithm  based 
on  de  Valk  [1990]  and  Rockafellar  [1974]  was  applied  to  the  2m  level  ADCP  (p,  u,  v) 
time- series  in  one-hour  increments  to  obtain  the  mean  wave  direction,  peak  wave 
direction,  peak  period,  and  spectral  width  of  the  incoming  waves.  The  Appendix  outlines 
the  procedure  for  obtaining  the  wave  energy  distribution  from  the  MISO  (p,  u,  v)  data. 
This  energy  distribution  was  calculated  over  a  frequency  (f)  range  of  0.03  </<  0.15  Hz, 
and  direction  (0  )  bins  of  five  degrees.  The  mean  wave  direction  (0„^„. )  was  calculated 
from  the  energy  distribution  function,  E(f,Q),  according  to 


j  jE(/,0)0de 

_  /-^ _ 

0.15  n 

j  jE(/,0)d0d/ 

0.03 -jr 


(27) 
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The  peak  direction  (0^)  was  calculated  by  summing  the  energy  into  direction  bins 
and  locating  the  maximum  energy  content.  Conversely,  the  peak  frequency  (fpeatd  was 
determined  by  summing  the  energy  distribution  into  frequency  bins  and  locating  the 
maximum  energy  content.  The  spectral  width  (IT)  was  determined  by 


W  =  2. 


M, 


(28) 


where  the  moments  Mo,  Mi,  and  M2  at  the  peak  frequency  and  across  all  directions  are 


(29) 

(30) 

=  (31) 


as  defined  by  the  directional  positions  relative  to  the  center  of  mass  (0com)- 

Figure  15  shows  the  energy  distribution  across  frequency  and  direction  for  one 
hour  of  an  energetic  winter  day.  The  spectral  energy  is  highly  focused  at  a  direction  of 
approximately  090°  and  a  /peak  of  0.0576  Hz  (Tpeak  =  IVs).  The  spectral  width  for  this 
hour  is  -13.1°.  Figure  16  contains  the  energy  distribution  of  a  less-energetic  summer 
day.  Note  that  the  energy  distribution  is  spread  across  all  directions  in  the  0.0667  Hz 
frequency  bin  (15s  period),  with  a  peak  direction  identified  at  90°.  Of  particular  interest 
is  the  apparent  reflection  of  energy  at  270°,  180°  opposite  of  the  peak  direction. 
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Oirectonal  Wave  Spetlra  for  OOOOZ  on  Yearday  200007S 
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Figure  15.  Wave  energy  for  a  high-energy  winter  yearday  (2000078)  according  to  frequency 
and  direction.  Note  the  highly  focused  energy  concentration  at  90°,  and  the  low 
spectral  width. 


To  ensure  that  reflected  energy  was  not  included  in  spectral  width  calculations,  spectral 
widths  were  only  calculated  from  0p-9O°  to  0p-i-9O°.  The  spectral  width  for  this  day  given 
these  constraints  was  a  large  70°. 


Infragravity  wave  energy  was  estimated  from  the  corrected  2m  height  velocity 
component  time-series  for  each  hour  by  low-pass  filtering  the  velocity  time-series  with  a 
fourth  order,  digital  Butterworth  filter  with  a  30- s  cutoff  to  remove  wave  velocities  and 
retain  only  infragravity  motions.  RMS  quantities  due  to  infragravity  motion  were 
calculated  to  characterize  the  level  of  infragravity  forcing  found  at  the  MISO  site. 


The  wave  orbital  radius  {dr)  was  estimated  from  the  time  integration  of  the  wave- 


frequency  energy  by  applying  the 


operator  in  the  frequency  domain  and  integrating 


across  the  wave  frequencies: 
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Oirectonal  Wave  Spectra  for  1 700Z  on  Yearday  2000281 


Figure  16.  Wave  energy  for  a  low-energy  (note  the  scale  change  in  energy  magnitude  from 
figure  15)  summer  yearday  (2000281)  according  to  frequency  and  direction.  Note 
the  spread  of  energy  across  all  directions  and  the  large  spectral  width.  Also  of 
interest  is  the  reflection  of  energy  180°  opposite  the  peak  wave  direction  of  90°  at 
270°. 
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where  (O  =  Inf ,  fy  is  the  nyquist  frequency,  and  Sxx(f)  is  the  autospectrum  of  each  hour- 

long  velocity  component  time- series.  A  30- s  high-pass  cutoff  has  been  selected  here. 
The  wave  orbital  diameter  0o)  was  calculated  from  the  magnitude  of  both  velocity 
components’  radii: 
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where  dr,u  and  dr,v  are  the  cross-shore  and  alongshore  wave  orbital  excursion  radii, 
respectively. 

Sediment  transport  is  governed  by  (among  other  factors)  the  bottom  stress. 
Bottom  stress  can  be  divided  into  two  components:  form  drag  and  skin  friction,  where 
skin  friction  is  primarily  responsible  for  sediment  motion.  Using  the  RMS  and  mean 
velocity  magnitudes,  the  bottom  stress  was  calculated  from  both  wave  (Xw)  and  current 
(Xc)  components: 


(34) 


where  p  is  the  density  of  seawater,  Urms  and  Vrms  are  the  wave  RMS  component 
quantities,  Um  and  Vm  are  the  mean  current  component  quantities,  and/w  and/c  are  wave 
and  current  friction  factors,  respectively.  The  wave  friction  factor  is  modeled  after 
Jonsson  (1966)  as 


/w  =exp 
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where  D50  is  the  median  grain  diameter  of  0.260  mm.  The  current  friction  factor  used  in 
equation  12  is  estimated  as  =  0.003  following  Amos  et  ah,  1988. 

The  Shields  parameter  was  utilized  to  determine  the  point  when  sediment  motion 
was  initiated  and  to  predict  sediment  transport.  It  represents  the  ratio  of  the  bed  shear 
stress  to  the  immersed  weight  per  unit  area,  and  is  estimated  by 
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(36) 


e'  = 


p{s-l)gD^o 


where  s  is  the  ratio  of  sediment  grain  density  to  water  density  (2.65  for  quartz  sand),  and 
g  is  the  gravity  acceleration. 


B.  SCANNING  ACOUSTIC  ALTIMETER 

I.  Data  Acquisition 

A  scanning  Acoustic  Altimeter  (SAA)  (Stanton,  2000),  developed  by  the  ocean 
turbulence  research  group  at  the  Naval  Postgraduate  School,  continuously  mapped  the 
sandy  bed  in  a  Im  alongshore  by  1.5m  cross-shore  area  immediately  offshore  from  the 
MISO  frame.  This  high-resolution  altimeter  consisted  of  a  transducer  head  mounted  on  a 
2- axis  stepper  motor  that  directed  the  altimeter  transducer  toward  the  ocean  bottom  in 
alongshore  (x)  and  cross- shore  (y)  directions  offset  by  12°  from  the  ADCP  orientation 
(Figure  6).  An  onboard  microprocessor  processed  the  acoustic  data,  found  the  bed  range, 
and  controlled  the  motors  as  the  SAA  scanned  in  1°  increments  over  a  ±35°  span  in  each 
axis,  allowing  the  changing  bumps  and  ripples  on  the  bed  to  be  measured  even  in 
optically  opaque,  turbid  water  conditions.  The  SAA  began  its  sampling  by  starting  at  a 
centerline  and  stepping  rapidly  back  and  forth  cross- shore  while  incrementing  the 
alongshore  angle  at  the  end  of  each  cross- shore  pass.  Upon  completion  of  this  half  of  the 
area  scan,  the  SAA  returned  to  the  centerline  and  progressed  out  to  the  left,  sampling 
each  point  with  20  pulses  of  1.3  MHz  sound  waves.  By  oversampling  spatially  and  using 
an  Objective  Analysis  technique,  a  gridded  z(x,y)  map  and  an  error  map  of  center 
weighted  averages  of  data  within  a  6  cm  diameter  Gaussian  field  at  each  grid-point  was 
created.  This  resulted  in  a  regular,  smoothed  morphology  map  with  a  1.2  cm  increment 
in  X  and  y  with  a  3  cm  half  power  bandwidth.  Ripple  wavelengths  of  0.03  to  1.2  m  were 
resolved  in  these  maps.  Figure  17  shows  typical  spatial  sampling  positions  of  the  raw 
bed  elevations  before  gridding  using  an  Objective  Analysis. 
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The  sampling  time  for  each  map  was  typically  around  15-20  minutes,  with  the 
sample  frequency  determined  by  the  acoustic  turbidity  and  near-bed  sediment  motion.  At 
any  given  time,  if  there  were  large  amounts  of  suspended  sediment  or  air  bubbles  in  the 
water  column  preventing  a  bed  detection,  the  SAA  would  retry  three  times  before 
incrementing  to  the  next  scan  position. 


Typtcai  X/Y  psatons  of  Scanned  Altimeter  Sampfej 


Alongshore  Distance  (m) 


Figure  17.  An  output  grid  of  x/y  positions  from  the  Scanning  Acoustic  Altimeter  (SAA) 
following  processing.  Note  the  ripples  on  the  left  and  right  side  of  the  map, 
observed  by  the  relative  closeness  of  the  x/y  points. 
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2.  Data  Processing 

In  order  to  best  characterize  the  bed’s  geometry  from  the  SAA  maps  and  its 
response  to  wave  and  current  forcing,  it  is  necessary  to  distinguish  relic  ripples  from 
active  cases.  Ripples  can  be  present  on  sandy  beds  over  a  wide  range  of  flow  conditions 
extending  from  minimum  shear  velocities  necessary  to  initiate  sediment  motion  up  to 
water  flows  intense  enough  to  produce  sheet  flow  sand  motion,  eliminating  any  ripple 
formations.  As  the  wave  and  current  forcing  begin  to  increase  over  a  planar  bed,  sand 
grains  suspended  in  the  water  column  by  mean  currents  begin  to  settle  and  align  with 
each  other  in  ridges  perpendicular  to  the  direction  of  oscillatory  flow.  These  ridges 
increase  in  spacing  and  become  more  asymmetric  as  vortices  form  on  the  ridges’  lee  side 
under  increased  flow.  The  previously  planar  bed  has  now  been  transformed  into  a  field 
of  long-crested,  asymmetric  ripples  typically  observed  during  periods  of  moderate  wave 
and  current  forcing.  As  strong  forcing  conditions  subside,  and  the  flow  conditions 
significantly  decrease,  the  ripples  cease  to  grow  or  migrate  and  become  relict.  Once  the 
wave  energy  increases  again,  small  wavelength  ripples  begin  to  form  on  the  long 
wavelength  relic  ripples.  To  understand  the  relationships  between  the  forcing  and  bed 
response  during  these  changing  forcing  conditions,  it  is  necessary  to  correctly  isolate  the 
geometry  of  the  actively  growing  ripples,  which  are  typically  atop  large  unchanging  relic 
ripples.  The  following  method  was  utilized  to  make  this  distinction  and  correctly  identify 
the  wavelength  and  height  of  actively  growing  ripples. 

Using  a  two-dimensional  FFT,  a  ripple  height  variance  (rj^)  spectrum  was 
calculated  for  each  OA  grid  map  in  wavenumber  component  (kx,  ky)  space.  These  maps 
were  averaged  into  one -hour  time  intervals  to  coincide  with  the  wave  and  current  forcing 
statistics  calculated  from  the  BADCP  data.  Gaps  not  longer  than  six  hours  in  the  mapped 
T]  time- series  were  interpolated  across  using  a  linear  interpolation  scheme. 

9  9 

To  reduce  the  effects  of  noise  in  the  p  variance  spectrum,  all  p  values  below  a 
determined  threshold  level  were  rejected  and  removed  from  the  spectral  bins.  A  noise 
threshold  level  of  4x10"  m  was  determined  from  p  variances  at  high  wavenumbers 
during  quiescent  conditions  under  the  assumption  that  white  noise  existed  across  the 
entire  time- series.  Figure  18  displays  the  p  time- series  used  in  the  selection  of  the  noise 
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level  threshold.  Nine  bins  of  high  wavenumber  (14.1-16.5  rh^)  values  collected 
during  32  calm,  low-energy  days  (768  hours)  contributed  to  this  time-series.  An  active 
ripple  limit  (above  which  ripples  may  legitimately  be  forming)  was  conservatively 
chosen  at  0.75x10“  m  .  Any  p  values  above  this  limit  did  not  contribute  to  the  noise 
threshold  calculations.  The  points  below  this  limit  were  averaged  to  determine  the  noise 
floor  for  the  entire  mapped  p  time- series. 
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Figure  18.  p  time- series  used  to  determine  the  noise  threshold  level  for  the  p  variance 

spectrum  derived  from  the  Scanning  Acoustic  Altimeter  (SAA)  objective  analysis 

9  -7  9 

maps.  Excluding  p  values  beyond  an  active  ripple  limit  of  0.75x10“  m  ,  nine 
bins  of  high  wavenumber  (14. 1- 16.5  m“  )  p  values  collected  during  32  calm, 
low-energy  days  were  averaged  to  determine  the  noise  floor  of  4x10“  . 
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Following  the  removal  of  all  p  values  below  the  noise  threshold,  each  (kx,  ky)  p 

value  was  normalized  by  its  bin’s  corresponding  mean  across  all  time  in  (kx,  ky,  t)  space: 
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2 

norm 


A  slew  rate  (5)  for  each  (kx,  ky)  combination  was  then  calculated  by  forming  a 
ratio  of  differenced  normalized  T]  values  to  time: 


dt 

A  fourth  order  low-pass  digital  Butterworth  filter  with  a  six  hour  cutoff  was 
applied  to  these  slew  rates  to  retain  the  six  hour  and  longer  changes  in  spectral  energy  in 
order  to  detect  significant  temporal  changes  as  the  ripple  field  evolved. 

Active  ripple  growth  periods  were  identified  whenever  the  slew  rates  exceeded  a 
threshold  level  of  3.25  m  /hour.  Figure  19  provides  an  illustrative  example  of  this 
threshold  selection  using  the  time-series  of  one  wavenumber  (K  =  3.00  rri*).  The  3.25 
threshold  level  (Figure  19(c))  identified  four  periods  of  ripple  growth  during  days  1-2,  9- 
11,  14-15,  and  27-28.  Decreasing  the  threshold  level  to  2.25  and  1.25  (Figure  19(a,b)) 
increased  the  number  of  active  ripple  detections  by  accepting  more  closely- spaced 
successive  r\l„^  values.  Conversely,  by  increasing  the  threshold  le’vel  to  4.25  (Figure 
19(d))  and  beyond,  large  differences  in  successive  values  were  required  for  active 
detection,  resulting  in  sparse  identifications  of  active  ripples.  For  these  reasons,  slew 
rates  exceeding  a  threshold  of  3.25  were  used  to  define  ripple  growth  for  each  ri^(kx,  ky). 

To  ascertain  the  ripple  geometry  during  these  periods  of  growth,  the  variance 
spectrum  was  binned  by  direction  into  alongshore  and  cross- shore  sectors  as  well  as  by  a 
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Aclivo  Ripple  Growlh  D«temfwiato«i  for  the  K  -  3.00  Normalsed  T«ne-Senee 
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Figure  19.  Threshold  selection  for  discerning  active  ripple  growth  from  relic  ripple  states  as 

9  - 1 

applied  to  the  normalized  p  variance  time-series  of  K  =  3.00  m'  . 
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radial- binned  wavenumber.  Mapped  values  within  the  alongshore  (0-45°)  and  cross¬ 
shore  (45°-90°)  sectors  were  then  binned  into  fifteen  evenly  spaced,  weight- averaged, 
composite,  radial  wavenumber  bins  ranging  from  0.83  m'^  to  11.67  rri^  Time- series  of 

'•y 

alongshore  and  cross- shore  p  variance  for  fifteen  wavenumbers  were  obtained  in  this 

'y 

manner.  Time- series  of  p  values  identified  as  active  according  to  the  filtered  slew  rates 

'y 

were  also  obtained.  If  an  p  value  was  identified  as  being  an  active  point  within  a 
direction  bin’s  wavenumber  bin,  all  p  values  in  that  wavenumber  bin  were  removed  and 
summed,  thus  producing  fifteen  p  sums  (0  in  place  of  inactive  cases).  From  this 
wavenumber  distribution  of  p  ,  the  ripple  wavelength  and  height  can  be  calculated 
according  to  the  weighted  centroid  of  the  maximum  variance  contribution  and  its  two 
adjacent  variance  values. 
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IV.  RESULTS  AND  ANALYSIS 


The  wave  and  current  forcing  time- series  for  year  2000  are  shown  in  Figure  20. 
Mean  currents  in  the  alongshore  direction  (not  shown)  with  peaks  of  up  to  22  cm/s  are 
stronger  than  the  weaker  tidal  currents  in  the  cross- shore  direction  (Figure  20(a))  with 
peaks  up  to  8  cm/s.  The  cross- shore  RMS  wave  velocities  (Figure  20(b))  reached  peaks 
of  20-30  cm/s  during  the  winter  (Days  1-100),  reducing  to  5-10  cm/s  during  the  spring 
and  summer  months  (Days  101-250).  Peak  periods  (Figure  20(d))  reached  13-20s  during 
the  stormy  winter  months,  and  10- 15s  during  calmer  summer  conditions.  The  spectral 
widths  (Figure  20(e))  of  incident  wave  energy  at  the  MISO  site  changed  from  narrow- 
band  ranges  averaging  20°  during  the  winter  to  35°  spreads  during  the  summer.  Wave 
orbital  velocities  (Figure  20(f))  reached  peaks  of  l-2m  during  periods  of  high  oscillatory 
flow  in  the  winter,  and  reduced  to  <  0.5m  with  decreased  wave  forcing  during  ihe 
summer.  The  Shields  parameter  representing  the  combined  wave  and  current  flow’s 
ability  to  move  the  sediment  and  induce  ripple  formation  is  plotted  in  Figure  20(g).  A 
critical  Shields  parameter  of  =0.04  (Amos  et  al.,  1988)  shows  an  approximate 
threshold  level  for  sediment  motion.  The  Shields  parameters  exceeded  0^^,^ during 
periods  of  intense  wave  and  current  forcing  in  the  winter  with  values  ranging  up  to  0.15. 
During  the  summer,  the  bottom  stress  due  to  waves  was  less  than  during  winter 
conditions,  producing  Shields  parameters  <  0.07  and  largely  below 

Using  these  forcing  characteristics  as  an  indication  for  the  possibility  of  ripple 
formation  and  upon  inspection  of  the  alongshore  and  cross-shore  active  wavenumber 
binned  variance  time- series,  active  ripple  conditions  were  identified  frequently 
through  the  winter  days,  while  relict  ripples  dominated  the  more  quiescent  summer 
periods.  The  focus  of  this  investigation  is  to  understand  the  processes  that  transform  the 
bed  by  actively  building  ripples.  For  this  reason,  three  representative  cases  of  evolving 
wave  and  current  forcing  events  from  the  winter  time  period  are  presented  and  analyzed. 
A  representative  example  from  the  summer  is  also  provided  for  contrast. 
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Figure  20.  MISO  time-series:  (a)  cross-shore  mean  currents,  (b)  cross-shore  wave  velocities, 
(c)  mean  wave  direction,  (d)  peak  period,  (e)  spectral  width,  (f)  wave  orbital 
diameter,  (g)  shields  parameter,  and  (h)  cross-shore  infragravity  RMS. 
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Wave  and  current  forcing  time- series  from  days  50-54  (Figure  21)  demonstrate 
the  potential  for  ripple  growth  under  increasing  wave  forcing  conditions.  Between  days 
50-52,  low  orbital  velocities,  high  wave  orbital  diameters,  and  sub-critical  Shields 
parameters  were  observed.  The  direction/wavenumber  binned  p  variance  time- series  at 
the  start  of  this  time-series  (50-50.3)  indicate  that  the  ripple  field  grew  to  be  dominated 
by  cross-shore  0.6  m  to  1.2  m  orbital  ripples  (Figure  22(a)).  Ripple  wavelength 
predictors  proposed  by  Nielsen  [1981]  (equations  1  and  5),  and  Clifton  and  Dingier 
[1984]  (equation  10)  illustrated  similar  bed  responses  to  the  forcing  conditions  during  this 
time  (Figure  29,  upper  panel).  Both  the  wave -dominated  Wiberg  and  Harris  [1994] 
model  (equations  10-15)  and  the  combined- flow  Khelifa  and  Oullet  [2000]  model 
(equation  16)  predicted  lower  wavelengths  below  0.4  m.  Between  50.3  and  52.4,  the 
wave  forcing  reduced  below  0'„,,  leaving  relict  ripples  with  these  wavelengths  and 
heights. 

Between  days  52.4  and  52.6,  there  was  rapid  alongshore  growth  in  0.2-0.6m 
ripple  wavelengths  as  the  orbital  diameter  rapidly  increased  to  the  sub -orbital  range. 
Interestingly,  the  0'  was  just  sub-critical  during  this  time  period,  but  clearly  the  bed  was 
active.  As  the  Shields  parameter  transitioned  above  critical  during  day  52.6,  active  ripple 
growth  ceased  in  the  alongshore  direction,  while  the  pre-existing  0.6  m  relict  orbital 
ripples  in  the  cross- shore  direction  were  decreased  in  height,  with  0.4  m  ripples  growing 
in  their  place.  This  transition  occurs  at  the  start  of  predicted  anorbital  motion  (Figure 
22(c)).  Khelifa  and  Oullet’ s  [2000]  combined  flow  model  provides  the  closest  prediction 
estimates  of  wavelengths  around  0.4  m,  while  Nielsen’s  [1981]  wave -dominated  models 
and  Wiberg  and  Harris’  [1994]  predictors  estimated  wavelengths  around  0.8  m. 

Days  28-31  provide  a  similar  example  of  alongshore  ripple  growth  (Figure  24(b)), 
but  with  a  more  gadual  and  lower  increase  in  wave  energy.  Throughout  the  first  day, 
there  is  sporadic  growth  in  both  the  alongshore  and  cross-shore  directions  at  0.2,  0.24, 
and  0.3  m  wavelengths  during  a  period  of  marginally  sub-critical  0' .  At  29.1-29.5,  there 
was  rapid  alongshore  growth  as  the  wave  forcing  gradually  increased  from  0.1  to  0.15 
m/s  (Figure  23)  and  the  wave  orbital  diameter  entered  the  suborbital  range  (Figure  24(c)). 
The  alongshore  growth  at  this  time  was  dominated  by  0.4  m  X  ripples  with  smaller  levels 
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of  growth  at  0.2  and  0.24  m  X.  Growth  around  0.4  m  wavelengths  was  predicted  by 
several  models  including  Khelifa  and  Outlet  [2000]  and  Clifton  and  Dingier  [1984] 
(Figure  30).  As  in  the  previous  example  of  days  50-54,  alongshore  growth  terminated 
upon  the  transition  above  around  day  29.6  of  the  28-31  time  period.  From  day  29.6- 
30.6,  remains  above  critical  while  cross- shore  0.6  and  1.2  m  ripple  growth  is 
observed,  coinciding  with  a  rapid  decrease  of  short  wavelength  ripples  in  the  alongshore 
sector.  Interestingly,  the  large  1.2  m  X  alongshore  ripple  energy  remained.  At  this  point 
in  the  time- series,  the  observed  wavelengths  depart  from  those  predicted  by  Khelifa  and 
Outlet’s  [2000]  model,  and  more  closely  resemble  those  predicted  by  Wiberg  and  Harris’ 

[1994]  suborbital  X  equation  (15)  (Figure  30,  upper  panel)  as  the  ratio  for  this 

/  ^50 

period  of  change  fluctuates  between  suborbital  and  anorbital  conditions. 

Days  75-77  provided  an  example  of  a  very  slow  wave  forcing  increase  where  0' 
remained  at  the  critical  level  (Figure  25).  Ripple  growth  was  observed  in  both 
alongshore  and  cross- shore  sectors  (Figure  26(a,b))  during  largely  orbital  conditions 
(Figure  25(c)).  Similar  to  the  previous  two  cases,  increases  in  alongshore  0.2-0.4  m 
ripples  were  observed  (days  75.4-75.7).  Khelifa  and  Outlet  [2000]  adequately  predict 
this  increase  under  combined  flow,  while  more  rapid  growth  is  estimated  by  Nielsen 
[1981]  and  Wiberg  and  Harris  [1994]  (Figure  31).  Unlike  days  50-54  and  28-31, 
however,  this  alongshore  growth  occurred  during  conditions  where  0 '  >  Qf, .  Cross¬ 
shore  0.6  and  1.2  m  ripple  growth  dominated  the  binned  spectrum  during  peak  orbital 
velocities  of  0.14  m/s  during  75.7-76.2,  however,  and  these  ripples  were  slowly 
decreased  between  76.2-88  as  orbital  velocities  and  displacements  lowered  slightly. 
Nielsen’s  [1981]  field  data  equation  (5)  and  Wiberg  and  Harris’  [1994]  suborbital 
equation  (15)  provided  the  closest  predictions  of  this  evolution.  Khelifa  and  Outlet’s 
[2000]  as  well  as  Clifton  and  Dingier’ s  [1984]  wavelength  equations  (16  and  8, 
respectively)  overestimated  the  bed’s  response  to  wave  and  current  forcing. 

Ripples  produced  during  relatively  calm  conditions  of  a  representative  summer 
period  provide  an  interesting  contrast  to  the  actively  changing  bedforms  indicative  of 
strong  winter  forcing.  Days  270-274  illustrate  a  long  period  of  ripples  in  a  relict  state  as 
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typically  observed  for  the  summer  months  (Figure  28).  As  the  orbital  motion  begins  to 
increase  and  0'  transitions  above  critical  (Figure  27),  alongshore  ripples  ranging  from 
0. 3-0.4  m  were  produced  around  day  274  and  275.4.  All  wavelength  predictors  initially 
provide  estimates  of  relic  ripple  states  during  days  270-272  (Figure  32).  At  day  272,  a 
slight  increase  in  wavelength  is  predicted  for  all  models.  During  the  period  of  active 
ripple  growth  (day  274-275.4),  Khelifa  and  Outlet  [2000]  once  again  predict  wavelengths 
comparable  to  those  observed  from  the  MISO  time- series.  The  Clifton  and  Dingier 
[1984]  and  Wiberg  and  Harris  [1994]  models  produce  similar  results,  while  Nielsen’s 
[1981]  ripple  wavelength  predictor  overestimates  wavelengths  by  0.2m. 
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Wave  and  Current  forcing  time-series  for  days  50-55:  (a)  Shields  paramet 
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Figure  22. 


Bed  response  time- series  for  days  50-55:  (a)  wavenumber  binned  r]  (cross¬ 
shore),  (b)  wavenumber  binned  (alongshore),  and  (c)  ripple  type  classification. 
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Figure  23. 


Yeiirchiy  -  2000 


Wave  and  Current  forcing  time- series  for  days  28-31:  (a)  Shields  parameter, 
(b)  wave  orbital  diameter,  (c)  wave  orbital  velocities,  (d)  mean  currents,  and 
(e)  infragravity  RMS  velocity. 
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Figure  24.  Bed  response  time- series  for  days  28-31 :  (a)  wavenumber  binned  r]  (cross¬ 
shore),  (b)  wavenumber  binned  (alongshore),  and  (c)  ripple  type  classification. 
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Figure  26.  Bed  response  time-series  for  days  75-77:  (a)  wavenumber  binned  r]  (cross¬ 
shore),  (b)  wavenumber  binned  (alongshore),  and  (c)  ripple  type  classification. 
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Wave  and  Current  Forcng  Time-Sehes  Summary 


Figure  27.  Wave  and  Current  forcing  time-series  for  days  270-276:  (a)  Shields  parameter, 
(b)  wave  orbital  diameter,  (c)  wave  orbital  velocities,  (d)  mean  currents,  and 
(e)  infragravity  RMS  velocity. 
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Figure  28. 
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Figure  29.  Bed  response  algorithms  for  (a)  ripple  wavelength,  and  (b)  ripple  height  according 
to  C&D  -  Clifton  and  Dingier  [1984],  N  -  Nielsen  [1981],  W&H  -  Wiberg  and 
Harris  [1994],  and  Khelifa  and  Gullet  [2000]  for  days  50-54. 
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Figure  30.  Bed  response  algorithms  for  (a)  ripple  wavelength,  and  (b)  ripple  height  according 
to  C&D  -  Clifton  and  Dingier  [1984],  N  -  Nielsen  [1981],  W&H  -  Wiberg  and 
Harris  [1994],  and  Khelifa  and  Gullet  [2000]  for  days  28-31. 
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Figure  31.  Bed  response  algorithms  for  (a)  ripple  wavelength,  and  (b)  ripple  height  according 
to  C&D  -  Clifton  and  Dingier  [1984],  N  -  Nielsen  [1981],  W&H  -  Wiberg  and 
Harris  [1994],  and  Khelifa  and  Gullet  [2000]  for  days  75-77. 
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Figure  32.  Bed  response  algorithms  for  (a)  ripple  wavelength,  and  (b)  ripple  height  according 
to  C&D  -  Clifton  and  Dingier  [1984],  N  -  Nielsen  [1981],  W&H  -  Wiberg  and 
Harris  [1994],  and  Khelifa  and  Gullet  [2000]  for  days  270-276. 
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V.  CONCLUSION 


The  analysis  of  bedform  evolution  under  combined  wave  and  current  forcing  at 
the  nearshore  MISO  site  has  detected  a  significant  amount  of  ripple  evolution  in  the 
alongshore  direction.  It  has  been  demonstrated  that  the  bed  is  dominated  by  relict  ripples 
remaining  from  previously  formed  ripples  as  the  forcing  decreased  below  critical  flow 
requirements.  As  the  combined  flow  increased  above  critical,  a  rapid  transition  from 
alongshore  to  cross- shore  bedform  evolution  was  observed  as  the  bedforms  adjusted  to 
new  active  forcing.  The  resulting  ripple  growth  was  not  wave  direction  oriented. 

Using  the  bed  geometry  and  wave  and  current  forcing  time-series,  it  has  been 
demonstrated  that  upon  comparison  of  several  laboratory  and  field  derived  models, 
Khelifa  and  Outlet  [2000]  provided  substantially  better  estimates  of  ripple  wavelength, 
even  during  periods  of  alongshore  growth.  Nielsen’s  [1981]  predictors  greatly 
overestimated  bed  response,  while  the  Wiberg  and  Harris’  [1994]  suborbital  model 
performed  adequately  during  periods  of  cross- shore  ripple  growth. 

To  better  understand  the  conditions  resulting  in  alongshore  ripple  growth  and 
account  for  this  growth  not  being  adequately  predicted  by  laboratory  models,  an  on-going 
investigation  will  examine  contributions  from  mean  currents,  infragravity  energy, 
reflected  wave  energy,  and  wave  asymmetries.  Statistical  analysis  of  these  time- series 
may  yield  important  insight  into  the  development  of  ripples  contrary  to  the  orientation  of 
wave  orbital  velocities. 
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APPENDIX 


To  calculate  the  directional  spectra  of  incoming  waves  at  the  MISO  site,  a 
Maximum  Entropy  Method  (MEM)  directional  wave  spectra  algorithm  was  applied  to 
pressure  ^),  cross-shore  ii)  and  alongshore  4')  velocity  time-series.  The  procedure 
utilized  is  based  on  de  Valk  [1990]  and  Rockafellar  [1974].  Eirst,  cross- spectral  densities 
are  obtained.  These  cross- spectral  densities  can  be  viewed  as 


4>i  ={pp) 

(Al) 

(|)2  =  (pu)  ~  ppcosQ 

(A2) 

^3  “  (t'^)  ~  PPsinO 

(A3) 

(|)4=(mm)~  ppcos^Q 

(A4) 

^5  =  ~  PP  sin  ^  6 

(A5) 

^6  =  ~  PP  sin  9  cos  6 

(A6) 

Using  sinusoidal  identities, 


<1^4  “<1^5  ~  pp|cos^0  -sin^0|  =  ppcos29 

(A7) 

2(|)g  ~  pp2sin0cos0  =  ppsin20 

(A8) 

Normalizing  the  Eourier  components  A1-A3  and  A7-A8  according  to  Long  [1980] 

yields 
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4*3  = 


4*3 


4*1  -3/4*2  “*“4*3 

/  _  4*4  ~4*5 


4*4  ~ 

(^5'  = 


4*4  +4*5 
2(^6 

4*4  +4*5 


(A9) 

(AlO) 

(All) 

(A12) 

(A13) 


To  calculate  the  entropy  of  system  (E(d))  as  a  function  of  the  directional 
distribution  (D)  on  (-7C,7l),  let 


n 

E{e)  =  -\D]n{D)dQ 

-K 


(A14) 


In  order  to  maximize  the  entropy  of  this  system  while  keeping  it  constrained  by 
the  measured  quantities  (normalized  spectral  values,  (|)'„,  calculated  in  A9-A13),  the 
equation 


n 

I d{B  )exp  {inB  )dB  -(|) '  =  0 

-71 


(A15) 


must  be  satisfied.  Introducing  the  Lagrangian  multipliers  p„  yields 
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Jl 

F=j  -D(e)]n(D(e))- J^^l„D(e)exp(me)  de  =0  (A16) 


Maximizing  F  with  respect  to  D  (i.e.  maximizing  the  entropy)  produces  the 
equation 


dF 

Id 


- 1  -  In  (D(e ))  -  |i„  exp  (me ) 


de  =0 


(A17) 


where 


In  (D(e ))  =  -1  -  |l„  exp  (me ) 


(A18) 


and  the  definition  of  directional  distribution  is 


D(e )  =  exp 


-1-Ll^«exp(me) 


(A19) 


Substituting  D  (A19)  into  F  (Alb), 


ti 

<7=  j-exp  -l-^|i„exp(me  )  *  -l-J^M-n  exp(me) 


+ 


(A20) 


-  1  -  E  )  ®^P  +  E 


and  simplifying 
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(A21) 


ji 

<7=  jexp  -l-^|i„exp(me ) 


To  maximize  q  as  a  function  of  the  Lagrangian  multipliers,  this  equation  is  solved 
by  the  BFGS  Quasi-Newton  method  with  a  mixed  quadratic  and  cubic  line  search 
procedure  (Broyden,  1970  and  Shanno,  1970). 

Applying  the  linear  dispersion  relationship  and  a  linear  theory  wave  transfer 
function 


H{k)  = 


^cosh{kh)^ 

cosh(^z)^ 


(A22) 


where  k  is  the  wavenumber  ),  h  is  the  local  water  depth,  and  z  is  the  depth  of 

measurement,  to  the  maximized  directional  distribution  multiplied  by  the  cross-spectral 
densities  produces  the  wave  directional  spectral  estimates  corrected  to  the  surface 
elevation.  From  this  quantity,  the  peak  frequency  (or  period),  peak  wave  direction,  mean 
wave  direction,  and  spectral  width  can  be  calculated. 
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